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ON  THE  CORRELATION  FUNCTIONS  IN  TIME  AND  SPACE 
OF  WIND-GENERATED  OCEAN  WAVES 


By 

Jan  Geert  de  Boer 

ABSTRACT 


When  wind-generated  ocean  surface  waves  are  described  statistically 
by  means  of  a  wave  energy  spectrum,  the  correlation  function  for 
the  surface  elevation  at  two  points  and  two  instants  of  time 
follows  as  an  integral  over  the  wave  spectrum.  The  time- 
correlation  function  and  the  spatial- correlation  function, 
which  are  important  for  the  statistical  description  of  an 
underwater  sound  field  scattered  from  the  sea  surface,  follow 
from  this  integral  as  special  cases.  They  are  examined  here 
for  two  proposed  spectra  for  fully-developed  seas  —  the  Neumann 
and  the  Pierson-Moskowit z  spectra  —  by  numerical  integration. 

It  is  shown  that  the  spatial- correlation  function,  although 
a  function  of  two  variables,  can  be  expressed  in  terms  of  two 
functions  of  only  one  variable  each,  when  a  cosine-squared 
law  for  the  directionality  of  the  wave  spectrum  is  assumed. 

These  functions  are  tabulated  and  plotted,  A  very  simple 
relation  is  sufficient  to  reconstruct  the  entire  anisotropic 
spatial  correlation  function  from  these  basic  functions. 
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INTRODUCTION 


The  subject  of  scattering  and  reflection  of  sound  waves  from 
rough  boundaries,  such  as  the  sea  surface,  has  received  increasing 
interest  in  the  past  fifteen  years.  Many  models,  both  of 
deterministic  and  random  character,  have  been  proposed  to 
describe  the  scattering  phenomenon.  A  detailed  discussion  of 
the  existing  literature  can  be  found  in  Ref.  1. 

The  elevation  and  slopes  of  the  sea  surface  are  random  processes 
in  time  and  space.  Consequently,  a  realistic  description  of 
the  scattered  sound  field  must  also  be  stochastic  in  nature. 
Statistical  quantities  such  as  mean  value,  covariance  and 
correlation  are  hence  of  interest;  and,  as  they  involve  the 
statistics  of  the  sea  surface,  knowledge  of  this  boundary 
is  required. 

If  the  statistical  description  of  the  scattered  field  is  limited 
to  first  and  second  moments,  the  sea  surface  is  sufficiently 
characterized  by  mean  value  and  covariance.  Moreover,  if  the 
mean  value  is  made  zero,  which  can  be  done  without  loss  of 
generality  (the  mean  level  does  not  change  in  the  time  interval 
that  is  typical  for  an  acoustical  experiment),  it  is  only  the 
correlation  function  of  the  surface  that  has  to  be  known. 

The  most  realistic  way  to  obtain  this  correlation  function  is 
by  using  the  (Neumann)  theory  of  the  surface  wave  spectrum  for 
a  fully-developed  sea.  This  spectrum  and  the  correlation  function 
are  related  via  the  familiar  Fourier  transform. 

Several  formulae  have  been  proposed  for  the  spectrum  of 
wind-generated  ocean  waves  in  a  fully-developed  sea.  It  is  the 
aim  of  this  report  to  compute  the  time-correlation  function  and 
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the  spatial-correlation  function  for  the  two  spectra  proposed  by 
Neumann  and  by  Pierson  and  Moskowitz.  It  should  be  noted 
that  the  Neumann  spectrum  has  now  been  discredited  on  both 
theoretical  and  experimental  grounds,  but  it  is  included 
here  for  comparison  and  because  it  is  so  commonly  referred  to 
in  the  existing  acoustic  literature  on  scattering  from  the 
sea  surface. 
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1. 


STATISTICS  OF  THE  OCEAN  SURFACE 


1 . 1  Introductory 

When  wind  is  blowing  over  the  surface  of  the  sea,  a  very 
complicated  mechanism  of  interaction  between  air  and  water  causes 
the  formation  of  surface  waves.  Many  studies  have  been  made  to 
investigate  this  phenomenon,  and  many  models  have  been  proposed 
to  describe  it.  But  a  description  that  covers  all  aspects 
is  not  yet  available. 

Attempts  have  been  made  to  characterize  the  sea  surface  with 
only  one  parameter,  especially  the  wind  speed.  But  the  time 
during  which  a  certain  constant  wind  has  been  blowing  (the 
"duration11)  and  the  size  of  the  area  over  which  it  has  been 
blowing  (the  "fetch11)  also  play  an  important  role.  This  has 
lead  to  the  concept  of  a  "fully-developed  sea",  over  which  the 
wind  speed  and  direction  have  been  constant  long  enough  for 
the  wave  system  to  contain  the  maximum  amount  of  energy 
it  can  possibly  have:  an  equilibrium  has  been  reached.  Clearly, 
this  is  only  a  theoretical  construction:  winds  of  constant 

speed  and  direction  do  not  last  very  long,  certainly  not  in 
large  areas.  Nevertheless,  the  idea  of  a  "fully-aroused  sea" 
has  produced  useful  results. 

A  very  good  introduction  to  the  subject  is  given  by  Kinsman 
[Ref.  21.  More  recent  insights  are  presented  by  Phillips  [Ref.  3]* 
Both  authors  point  out  that  the  sea  surface  is  a  random  process, 
in  space  as  well  as  time.  This  process,  z  =  ^(x,  t),  is  not 
Gaussian  (there  is  a  certain  skewness  of  the  waves,  and  waves 
of  infinite  height  have  zero  probability),  but  in  many  respects 
it  may  be  considered  as  "quasi-Gaussian " ,  as  measurements  have 
indicated. 
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In  its  most  general  form  the  second  moment  of  the  process  is 
the  correlation  function 

F  “  E  [  C(x^  ,  t^)  C(x2  >  t2)]  *  [Eq.  1*] 

In  principle  it  depends  on  the  position  of  the  points  of 
observation  and  on  the  observation  times.  But  the  assumption 
that  the  sea  surface  is  homogeneous  (at  least  in  the  area  where 
an  experiment  takes  place)  and  stationary  (at  least  for  the 
duration  of  an  experiment)  reduces  the  correlation  function  to 
a  function  of  only  space  differences  and  time  differences: 

E[£(x  t  )  C(x  ,t,  )]  sE[£(x  ,t  )  G(x  +  p,  t  +  T)] 

11  C  d  ll  1  1 

=  F(p>  t)  , 

where 

-4  - 4 

p  =  x .  -  x^  and  t  =  tg  -  t^  . 

In  the  past  it  has  sometimes  been  assumed  that  the  spatial- 
correlation  function  of  the  sea  surface  had  the  shape  of  a 
Gaussian  curve,  but  this  assumption  has  turned  out  to  be  incorrect 
and  unsatisfactorily.  r,The  most  realistic  way  to  incorporate 
the  correlation  function  of  surface  height  and  slopes  is  via 
the  theory  of  the  surface  wave  spectrum11  [Ref.  1,  p.91]* 

According  to  this  theory  the  sea  surface  is  considered  as  the 
combined  effect  of  a  large  band  of  sinusoidal  surface  waves  that 
travel  over  the  surface  in  very  many  directions,  each  having  its 
own  speed  and  hence  its  own  wave  number.  In  this  way  the  idea 
of  a  surface  wave  energy  spectrum  has  been  formed. 

Many  waves  travel  in  the  direction  of  the  mean  wind.  But  waves 
are  also  generated  sideways,  at  least  in  the  angular  interval 
(-tt/2,  rr/2),  and  may  be  even  in  an  interval  (-rr,  rr)  ,  because  of 

non-linear  wave  interactions.  It  can  be  said,  however,  that 
waves  are  generally  strongest  in  the  down-wind  direction.  Hence 
the  sea  surface  is  anisotropic,  and  the  wave  spectrum  depends 


[Eq.  2] 

[Eq.  3  ] 
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not  only  on  frequency  or  wave  number,  but  also  on  direction. 


1 . 2  Relations  Between  Surface  Correlation  Functions 
and  Wave  Spectra 

In  its  most  general  form  the  surface  correlation  function  deals 
with  the  surface  elevation  at  two  different  points  and  at  two 
different  instants  of  time.  Other  correlation  functions  can  be 
derived  from  it  by  taking  the  time  difference  equal  to  zero, 
by  considering  coinciding  points,  or  by  assuming  the  surface 
to  be  isotropic.  For  each  correlation  function  obtained  in 
this  way  a  corresponding  spectral  function  can  be  defined  via 
a  Fourier  transform  relation.  These  relations  are  discussed 
in  detail  in  Refs.  2  and  3* 


Returning  to  Eq.  2,  we  normalize  F: 

F(p,  T)  =  h2N  (p,  T),  [Eq.  4] 

h2  -  F ( 0 , 0 )  ;  [Eq.  5] 

the  vector  p  has  components  %  and  r\,  in  the  X  and  X 
directions  respectively. 


The  correlation  function  N("p,  t)  is  the  three-dimensional 
Fourier  transform  of  the  most  general  wave  spectrum  Y(l?,  a)  > 
a  function  of  both  wave  number  and  frequency: 


h2N (  p,  t) 


_  r 


1  ■+ 
dK 


w 

°  -4  -4-4 

d cf  Y(i?,  a1)  exp[i  (K*  p  -  c!  t)  ] 


[Eq.  6] 


For  small-amplitude,  deep-water  waves  there  is  a  unique  relation 
between  the  wave  number  and  the  frequency*  the  dispersion 
relation 


K  =  a2/ g  . 


[Eq.  7] 
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This  implies  [Refs.  2  and  3]  that 


Y(K,  a'  )  =  Y(K)  6(0  -  O'  )  ,  [Eq.  8] 

so  that 

f p  -♦  -♦  -4  -> 

h2N  (  p,  t)  =  jj  dK  Y(K)  exp[K  -p  -  ot)] 

—  00 

follows.  With  polar  coordinates  (K,  6)  and  again  the  dispersion 
relation,  we  then  find 


[Eq-  9] 


TT 


h2N( n)  = 


d0 


TT 


da 

o  g 


2  a 


[y(k,  e)] 

K=a2/ g 


Putting 


exp 


i  -2—  (  ^  cos  9  +  ri  sin  0)-iaT 

-  B  — 


[Eq.  10] 


*(a,  0)  =  —  [Y(K,  0)] 

g2  K=a2/g 


[Eq.  11] 


into  Eq.  10,  and  realizing  that  N  is  a  real  function,  we  get 
f inally : 


N(  ?,  n,  t)  =  (  2h 2  ) 


-1 


TT 


P 

J_  TT 


de 


da  $(  a,  0)  cos 


- —  (  5  cos  0  +  n  sf  n  8) 

g 


[Eq.  12] 


This  formula  is  the  starting  point  for  our  calculations.  It  may 
be  compared  with  Ref.  2[p„378,  Eq. (8.3:3)]- 


The  spectral  function  $(  a,  0)  is  non-negative  over  the  entire 
(  <j,  0)  -  plane.  The  integral 


I 


02  a2 

i  d  9  ,  da  §(  c,  0) 

Jel  Jol 


[Eq.  13] 
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gives  a  measure  of  the  energy  in  the  wave  field  having  frequencies 
between  and  gg ,  and  travelling  in  directions  between  0.^ 

and  02 •  If  we  integrate  the  spectrum  over  the  whole  ( g,  0) -plane, 
we  have  an  estimate  of  the  total  energy  in  the  wave  field: 


E 


,  TT 


d  0 


L-rr 


P 

\  da 

J 

o 


$(  a,  e)  . 


[Eq.  14] 


1 . 3  The  Sea-Surface  Roughness  Spectra 


Most  oceanographic  literature  about  the  surface  wave  spectra 
deals  with  the  function  $(  a)  ,  rathern  than  with  §( a,  0) • 

They  are  related  as  follows: 

[Eq.  15] 


§(  a) 


TT 


d  0  §(  a,  e) 


-TT 


When  §(  a)  is  given  and  §(  c,  0)  is  needed,  an  inverse  relation 
is  required. 


There  is  some  disagreement  in  the  literature  about  the  explicit 
form  of  the  function  $( a) .  Part  of  the  discrepancies  can  be 
explained  if  we  realize  that  the  measurements  on  which  the 
empirical  formulae  for  $ ( CJ )  are  based  have  not  all  been  made 
in  seas  with  the  same  state  of  development  [Ref.  1,  p.88-89] • 

A  formerly-used  estimate  for  the  functional  form  of  $( a)  is: 

He)  =  \  C  a~(>  exp(-2gs  a-2  U-3  )  ,  [Eq.  16] 

where  C  =  3*05  m2/s~*.  This  is  the  Neumann  spectrum,  usually 
called  A2(q)  [see  Ref.  2,  p.  389]*  Its  anisotropic  version  is 
given  by  Pierson,  Neumann  and  James  [see  Ref.  2,  p*399~j: 

$( o,  e)  h a2( a,  e)  =  c  o_t1exP(-2g2  o"2u"F)  cos2e 

for  o,  s  a  <  «°,  |  6  |  ^  —  , 

=  0,  otherwise.  [Eq.  17] 
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Comparison  between  Eqs.  16  and  17  yields  a  relation  between  the 
anisotropic  spectrum  and  its  reduced  form,  which  will  also  be 
used  for  other  spectra: 


H  o,  e) 


COS2  0 


=  0,  otherwise  \ 


^  a  < 


,  |  e  |  <;  TT/  2 


[Eq.  18] 


Another  spectral  function  is  mentioned  by  Schulkin  [Ref,  4]; 
it  has  a  sharp  low-frequency  cut  off: 


Y(K, 

e>  =| 

Y(K)  cos2  0  for 

CD 

IA 

N>|=t 

=  0 

for 

1*1  >1 

[Eq. 

19] 

with 

r.  2TT 

Y(K) 

'I 

d  0  Y(K,  0) 

[Eq. 

20] 

-  4 

„  it  -2 

=  BK 

for  K  ^gU 

1  1 

it  —  2 

[Eq. 

21] 

=  0 

for  | K |  <  g  U 

5 

this  is 

the  Burling-Phillips  spectrum. 

The  dimensionless  constant  B 

equals  0.46  X  10 

according  to 

_ 2 

Schulkin,  and  0,6  x 10  according  to  Phillips  [Ref.  3] • 

The  expression  for  $( a) ,  as  given  by  Neumann  [Eq,  16],  is  most 
often  met  in  the  literature.  But  its  asymptotic  behaviour  at 
high  frequencies,  showing  a  proportionality  with  c  ^ ,  has  been 


Cox  and  Munk,  and  Wong,  have  found  from  ratios  of  up-wind  to  cross-wind 
mean  square  slopes  that  the  cosine-squared  law  is  too  narrow  to  describe 
the  directionality  of  the  wind  waves  (see  Ref„  4,  pc17)0  We  note  in  passing 
that  the  quantity  "mean- square  slope'1  is  an  important  measure  of  sea  surface 
roughness;  it  includes  the  roughness  contribution  of  the  small  wave  patches 
riding  atop  the  large  sea  waves,  which  are  important  for  acoustics  and 
radar o 
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proved  wrong.  Pierson  and  Moskowitz  [Ref.  5]  gave  arguments  in 
favour  of  a  a  ^  behaviour.  They  stated  that  "within  the  present 
limitation  of  data,  the  spectrums  of  fully  developed  wind-generated 
seas  for  winds  measured  at  19*5  metres  are  given  very  nearly  by" 
[Ref.  5,  p.5190]: 


$(  c)  =  ag2  <7  ^  exp 

with  a  =  8.10  X 10  3,  3  =0.74  and  a  =  g/U,  where  U  is 

o 

the  wind  speed  reported  by  the  weather  ships. 


■&y 


[Eq.  22] 


A  formula  in  which  the  wind  speed  does  not  appear,  and  that  also 
has  a  high-frequency  behaviour  of  the  type  a  has  been  suggested 
by  Phillips  [Ref.  3].  He  gave  the  formula 

$(a)  =  c»  g2  a-5  ,  [Eq.  23] 

t  —  2 

where  C  =  1.2x10  ,  a  dimensionless  constant,  that  varies 

slightly  with  fetch. 


In  this  report  we  will  concentrate  on  two  spectra,  both  proposed 
for  a  fully-aroused  sea,  namely  the  spectra  given  in  Eq.  16 
(Neumann)  and  Eq.  22  (Pierson-Moskowitz) .  Since  the  anisotropic 
version  is  required,  we  will  also  apply  Eq.  18,  with  0^=0. 


The  frequencies  a^,  at  which  the  energy  spectra  reach  their 
maximum,  and  the  maxima  themselves,  can  be  found  easily.  We  have: 


a )  Neumann  Spectrum 

CTM  =  U  =  °'816  u  ’ 

*(aM)  =  f  U6  ; 


[Eq.  24] 
[Eq.  25] 
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b)  Pierson-Moskowitz  Spectrum 

°M  “  V  =  °-8”  f  ^ 

f(  o^)  =  f  U5  .  [Eq.  27] 


We  see  that  with  increasing  wind  speed  decreases  and  $(  C3^) 

increases.  For  the  surface  wave  spectrum  this  means  more  low 
frequencies  and  so  a  larger  correlation  distance. 


A  fully-developed  sea  for  a  fixed  wind  speed  U  is  one  whose 
spectrum  contains  components  of  all  frequencies  0  ^  Q  <  »  ^ 
each  with  the  maximum  energy  of  which  it  is  capable  under  the 
given  wind.  The  total  energy  in  such  a  sea  is  [cf  Eqs.  14  and  15] 


E 


da  $(  a)  =  2h2  . 


[Eq.  28] 


So  we  get  for  the  spectra  under  study: 


a)  Neumann  Spectrum 


b )  Pierson-Moskowitz  Spectrum 


[Eq.  29] 


[Eq.  30] 


Details  of  the  calculation  can  be  found  in  Appendix  A. 
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2  „ 


THE  SURFACE  CORRELATION  FUNCTIONS  IN  TIME  AND  SPACE 


The  foregoing  discussion  of  the  surface  wave  spectra,  and  the 
relation  between  the  directional  spectrum  and  its  reduced  form, 
enable  us  to  rewrite  the  formula  for  the  surface  correlation 
function  N  ( "p,  t)  .  From  Eqs.  12  and  18,  with  =  0,  we  obtain: 

[Eq.  31] 

The  normalized  time-correlation  function  N(0,t),  and  the 
normalized  spatial-correlation  function  N("p,  0)  follow  immediately 
from  this  expression.  They  will  be  calculated  in  the  following 
sections . 


h2N(  p,  T)=  i  J 


tt/2  pp 

d0  ,  da  $(  a)  cos2  0  cos 
-tt/ 2  J 


(  ^cos  0  +  ns  in  0)  -  c 


For  computational  reasons  we  bring  Eq.  31  into  a  more  convenient 
form:  we  divide  by  h2,  change  the  order  of  integration,  and 

separate  the  variables  in  the  argument  of  the  second  cosine 
function.  Then  we  get: 


N(p,  t)=(  th2 ) 


with 

Ij(  a,  p)  = 


2  \  -  1 


—  0 


da  $(  a)  cos(ai)l1(a,  p)+  j  da  $(  c)  sin(  ot)  I2  (  a,  p) 

[Eq.  32] 


[Eq.  33a] 


O  tt/2 

2 

d 0  cos2  0  cos 

— —  (  £cos0  +  rising) 

0 

_g  _ 

-tt/2 


pn/2 

I  (  a,  p)=  d  0  cos2  9  sin 


J  . 

-tt/2 


—  (  ?cos  0  +  n  sin  0) 

s 


[Eq.  33b] 


We  may  note  in  passing  that  h2  can  be  calculated  with  Eq.  28: 


h2 


1 

2 


oo 

da  §(  a) 
o 


[Eq.  34] 
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an  expression  that  follows  also  from  Eq.  31  by  putting 
and  t  =  0 .  The  quantities  E  for  the  Neumann  spectrum  and 
the  Pierson-Moskowitz  spectrum  are  given  in  Eqs  .  29  and  30, 

respectively.  Hence  we  find  for  h2: 

a )  Neumann  Spectrum 


[Eq.  35] 


with  C  =  3.05; 


b )  Pierson-Moskowitz  Spectrum 


[Eq.  30] 


with  a  =  8 . 10  x 10  3 ,  3  =  0.74. 
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3  - 


THE  TIME -CORRELATION  FUNCTION 


3 . 1  Introductory 

The  correlation  function  N(p, T)  given  in  Eq.  32  is  too  general 
for  a  detailed  analysis.  In  this  section  we  want  to  study  its 
behaviour  as  a  time-correlation  function.  Experimentally,  this 
function  could  be  obtained  by  observing  the  wave  pattern  at 
only  one  point,  but  for  a  long  time.  The  time  record  from  a 
shipborne  wave  height  meter,  or  from  a  wavepole  [both  giving  an 
estimate  of  £(t)  ]  could  serve  as  the  raw  data.  The 
function  N(0,t)  could  then  be  calculated  approximately  with 
the  formula 


t0+T 

dt  C(t)  C(t +  t)  .  [Eq. 37] 

0 

From  the  theoretical  side,  as  is  our  concern  in  this  report, 

N(0,t)  for  a  fully-developed  sea  follows  from  Eqs.  32  and  33  by 
setting  p  =  0,  i.e.  ^  —  T|  =  0.  We  then  have 


N(0,t)=| 


N  ( 0,  t)  =  (  nh 2 ) 


z\~l 


da  $(  a)  cos(ot)  ( a, 0)  , 


with 


PU/2  2 

I]L  (  a,  0)  =  d9  cos  8. 

-tt/2 


[Eq.  38] 


[Eq.  39] 


But  the  integral  I ^  is  easily  evaluated;  its  value  is  Tr/2. 
And  so  we  have  finally 


N(0,  t)  =  (  2h2 )  1  J 

0 

which  says  that  N(0,t) 


da  $(  a)  cos(  0T)  , 
is  the  Fourier 


cosine 


transform  of 


[Eq.  40] 

$(  a) . 
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3 . 2  The  Neumann  Spectrum 


In  Eq.  40  we  have  to  substitute  the  energy  spectrum  given  by 
Eq,  16,  and  the  value  for  hz  from  Eq.  35*  The  final  result  can 
be  written  in  a  more  convenient  form  if  we  introduce  the 
variable 

TT  ^ 

z  =  ~  ,  [Eq.  41] 

2g 

and  normalize  the  time  delay  T  by: 

Tn  =  T/U.  [Eq.  42] 

We  note  in  passing  that  a  similar  normalization  can  be 
found  in  Ref.  5* 

With  these  modifications  we  get: 

N(0,  t„)  =— -  }  dz  z“32  exp(-i)  cos  (gj2  Jz  tn)  .  [Eq.  43] 

J0  z 

Clearly  the  wind  speed  U  does  not  influence  the  shape  of  the 
curve,  but  only  the  time  scale. 


Unfortunately,  we  have  not  been  able  to  solve  this  integral 
analytically.  Hence,  a  numerical  integration  has  been  performed. 
The  result  is  plotted  in  Fig.  1. 


FIG.  1  TIME  CORRELATION  FUNCTION  OF  WIND  GENERATED  OCEAN  SURFACE  WAVES 

The  normalized  t  me  equals  U~*  t  mes  the  actual  t:me 
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It  should  be  noted  that,  although  N(0,t)  and  §( a)  are  related 
via  the  familiar  Fourier  cosine  transform  [Eq.  40],  common 
numerical  techniques  such  as  the  fast  Fourier  transform  (FFT) 
cannot  be  used.  This  is  because  the  FFT  requires  an  equidistant 
sampling  of  $(c?),  which  is  not  suitable  for  this  function  since 
it  has  a  steep  rising  part  (0  ^  0*  ^  o^)  and  a  very  slowly 
descending  part  (  0^  ^  C  <  »)  , 


3 . 3  The  Pierson-Moskowitz  Spectrum 

In  this  case  the  spectrum  of  Eq.  22  and  h2  of  Eq.  36  have  to  be 
substituted  into  Eq.  40#  This  time  we  change  from  q  to  the 
variable 

z  =  e( a0/ a) 4 ;  [Eq.  44] 

the  result  is 

0 00  l  _jl 

N(0,  Tn)  =  j  dz  exp(-z)  cos  (  04  g  TN  z~4  )  ,  [Eq.  45] 

again  with  =  t/U.  This  integral  has  also  been  approximated 

numerically.  The  result  is  presented  in  Fig.  1,  together  with 
its  deviation  from  the  Neumann  curve. 


3 • 4  Discussion 


In  Ref.  6,  Latta  and  Bailie  used  very  sophisticated  mathematical 
techniques  to  calculate  analytically  the  Fourier  cosine  transform 
of  the  energy  spectrum  [see  Eq.  40]  for  the  Neumann  and  the 
Pierson-Moskowitz  spectra.  The  resulting  formulae  are  so 
complicated,  however,  that  they  also  had  to  use  a  computer  for 
the  final  evaluation  of  the  function  N(0,t).  Figure  2  shows 
a  copy  of  their  curves,  for  comparison  with  Fig.  1;  it  also 
includes  a  curve  found  by  Bendat  [Ref.  7] >  for  the  Pierson- 
Moskowitz  spectrum. 
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S  e  10  12  H  16  li  20 


FIG.  2  TIME  AUTO-CORRELATION  FUNCTIONS  ON  A  NORMALIZED  SCALE  ( F-om  Ref  6) 

Unfortunately,  Latta  and  Bailie  plotted  their  results  on  a 
normalized  time  scale,  without  giving  explicitly  the  normalization 
factor,  so  that  a  quantitative  comparison  between  their  curves 
and  ours  becomes  quite  difficult. 

The  first  thing  we  notice,  when  comparing  Figs.  1  and  2,  is  that 
if  we  multiply  our  time  scale  by  ten,  our  PM-curve  is  in  good 
agreement  with  that  of  Bendat,  who  used  a  very  simple  integrating 
routine.  Because  the  Latta-Bailie  curve  is  only  slightly 
different  (their  result  is  probably  more  accurate,  since  their 
calculation  is  far  more  rigorous  than  Bendat' s  and  ours),  we 
might  conclude  that  they  have  normalized  the  time  scale  in 
accordance  with  Ref.  5  [p. 51821,  i.e.  using  =  gT/U. 

However  this  would  include  a  normalized  frequency  U)  =  ~~  and 
Eq.  22  would  then  change  into  (in  the  notation  of  Ref.  6) 


[Eq.  46] 


with 


B  =  0  =  0.74 


[Eq.  47] 


Nevertheless,  Fig.  2  gives  the  value  B  —  1.25.  Hence,  — gT/U 

or  t„  =  t/yU  with  y  =  0.10  20,  does  not  seem  to  be  the 
N 
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normalization  formula  used,  but  it  is  very  close  to  it.  Indeed, 
if  we  solve  the  equation 

B(^)  =  Bll)  [Eq.  48] 

we  obtain 

0)  =  (§)*  Y  ,  [Eq.  49] 

and  with  B  =  1.25,  g  =  0.74  and  g  =  9*81  this  gives 

uo  =  YUo  ,  A  =  y  |  —  ,  [Eq.  50] 

g2 

with  y  =  O.II65.  Figures  1  and  2  can  hence  be  compared  if 
we  multiply  our  scale  by  y  1  =8.6.  A  difference  between  the 
curves  from  the  Neumann  spectrum  is  then  apparent,  and  this 
discrepancy  is  hard  to  explain,  because  of  lack  of  information. 

First  we  note  that  Latta  and  Bailie  have  also  normalized  the 
Neumann-spectrum  [Eq.  16] ;  they  wrote 

§(  ul))  =  Kuf6  exp(-AU)2)  .  [Eq.  51] 


The  parameter  A  has  not  been  defined.  With  UU  =  yU  a  we  would 
get 

JL 

A  =  2 y2  g2  =  2(|)2  ,  [Eq.  52] 

which  is  not  in  agreement  with  Eq.  50.  This  leads  us  to  the 
conclusion  that  the  symbol  A  has  been  used  for  two  completely 
different  quantities . 


In  the  final  formula  for  N(0,t),  calculated  for  the  Neumann-spectrum 
[Ref.  6],  the  parameter  A  appears  again,  as  it  should; 


n 


N ( 0 ,  t)  =~  J~n  £ 

J  n=0 


(-l)n(nzl)(H|3)  a2  Tn 

r(n+i)  r  (4-1) 


[Eq.  53] 
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Nevertheless,  in  Fig.  2a  curve  is  plotted  without  a  value 
assigned  to  A. 

For  all  these  reasons  we  shall  not  attempt  a  quantitative  comparison 
between  Figs.  1  and  2,  as  far  as  the  Neumann  curve  is  concerned. 

We  only  conclude  that  the  curves  show  a  similar  behaviour. 

The  fact  that  the  time  scale  in  Fig.  1  could  be  normalized  by 
dividing  the  real  time  delay  by  the  wind  speed  is  very  interesting. 
It  means  that  a  fully-developed  sea  does  not  change  its  shape,  but 
only  its  scale,  when  the  wind  speed  is  changed  and  a  new 
equilibrium  has  been  reached.  The  whole  wave  pattern,  in  a 
statistical  sense,  is  stretched  when  the  wind  speed  increases 
and  contracted  when  it  decreases.  A  fetch  of  infinite 

length  and  a  constant  wind  of  infinite  duration  are  essential 
for  this  interpretation. 
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4- 


THE  SPACE-CORRELATION  FUNCTION 


4  *  1  Introductory 

When  sound  waves  are  scattered  from  the  sea  surface  not  one  point, 
but  a  whole  surface  area,  is  involved.  This  explains  why  the 
spatial- correlation  function  N(p,  0)  is  important  for  the 
description  of  the  scattered  sound  field.  In  fact,  its 
importance  is  far  greater  than  that  of  N(0,t). 

Experimentally  N(*p,  0)  is  hard  to  obtain.  In  principle,  its 
calculation  requires  the  surface  elevation  in  a  large  area,  at 
one  instant  of  time,  i„e,  £(x,  t)  for  t  =  t0  and  x  ranging 
through  a  certain  area  A  large  enough  to  include  the 
biggest  p  of  interest.  Given  this  information,  N(*p,  0)  could 
then  be  computed  with  the  formula 

N(p,0)  =j  |  |  dx  C(x,t0)  C(x+  p,  t0)  .  [Eq.  54] 

An  approximation  made  using  a  set  of  sample  positions  5^.  .  is 

1 9  3 

also  possible: 


^  C(3i,j '  *«>  m>  *»>  ’ 


[Eq.  55] 


where 


l,rn 


=  x 


i+l,  J+m 


-♦ 

-  X  . 


J 


This  would  require  FQ  wave-height 


meters,  whose  positions  are  fixed  and  well  known,  and  a  synchronous 
recording  of  their  readings.  Such  an  experiment  is  hard  to 
imagine . 


More  promising  seems  the  optical  method  described  by  Stilwell 
[Ref.  8  j,  in  which  the  sea  surface  is  illuminated  by  a  "continuous 
skylight”  and  recorded  photographically.  Such  a  record  indeed 


21 


represents  a  continuous  surface  area  at  one  instant  of  time.  The 
data  processing  is  quite  simple.  !lBy  the  use  of  optical  analysis 
it  is  possible  to  resolve  the  variations  of  density  in  the 
surface.  When  a  transparency  of  a  surface  photograph  is  placed 
in  one  focal  plane  of  a  lens,  the  Fourier  transform  of  the 
variations  appear  as  light  amplitude  in  the  other  focal  plane. 

This  information,  in  addition  to  the  height  and  aspect  angle  of  the 
camera,  allows  the  energy  spectrum  of  the  surface  to  be 
obtained”.  [Ref.  8,  Abstract]. 


When  the  surface  wave  spectrum  is  given,  as  in  our  case,  N(p,0) 
can  be  calculated  from  Eqs.  32  and  33  by  taking  t  =  0 .  Doing 
so  we  get 


N(p,0)  =  (  Tth 2  ) 
with 


-1  n  00 


da  $(  a)  i^a,  p) 


I  x(a,3)  = 


TT/2 


-Tt/2 


dQ  cosfc Scos 


(  5cos  8  +  T)sin  0) 


[Eq.  56] 


[Eq.  57] 


The  integral  1a  can  be  expressed  in  terms  of  Bessel  functions 
[see  Appendix  B  for  details] : 


Tl  (  o,  p)=  rrp 


-2 


-1 


J 


>T 


)  +  (  n2  - 


5’X'T' 


V 


2 

-G- 

g 


[Eq.  58] 


4 . 2  The  Neumann  Spectrum 


In  Eq.  56  we  substitute  Eqs.  58,  16  and  35*  Moreover,  we  introduce 
a  normalized  correlation  distance  ,  with  components 

and  n^T,  by  using  the  relation  [cf  Ref.  5,  p.5l8  2]: 


N 


[Eq-  59] 


u* 


and  change  the  variable  a  by  taking,  as  in  Sect.  3.2: 


CgU  2 

2g2 


[Eq.  60] 
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Then  the  correlation  function  N  ("p,  0)  can  be  written  in 
Cartesian  coordinates  (with  p^T  =  )  *  as: 


N(V  ^  »  0)  =  P; 


N  |jN  K  1  'DN 


>  +  <V-Sn8>Mpn> 


[Eq.  61] 


or  in  polar  coor  dinates  (  cos  cp,  'Hjj  =  Pj,j  s;*-n  Cp)  as 


N(  PN,  Cp,  0)  -  cos2  CpKt  (  pjj)  +  ( sin 2  cp  -  cos2cp)  Ks  (  pN)  , 
where 

8 


K1<Pn>  = 


3,rs  Jc 


r*  dz  e-1/^  .-7/2 


J  o  (  PN  7  ) 


[Eq.  62] 


[Eq,  63a] 


Mpn>  " 


3  J~v 


P®  .  -1/z  -7/2  J1^PnZ-) 

,  dz  e  '  z  '  — 7 - 

J  (pNz) 


[Eq.  63b] 


The  functions  and  K2  have  been  calculated  numerically,  on  a 

digital  computer  (Elliott  503)*  and  are  plotted  in  Fig„  3*  For 
later  use  their  sample  values  are  given  in  Table  1.  The  fact 
that  they  are  one-dimensional  simplifies  the  calculation  of  the 
two-dimensional  function  N(f^,  *  0)  considerably:  we  only 

need  to  read  the  values  of  K1  and  Kp  into  the  computer, 
after  which  Eqs.  61  or  62  enables  us  to  compute  any  value  of 
the  surface  N(^,  t]^,  0).  Moreover,  N  is  symmetric  with 

respect  to  the  diagonal  =  'Hjj  : 


N(?n,  V  0)  =  N(V  V  0)>  CE(i*  64^ 

this  halves  the  number  of  samples  to  be  computed,. 

In  Fig,  4  we  have  plotted  the  correlation  function  N (  ^ ,  T)^ ,  0) 

for  the  Neumann  spectrum.  The  scales  are  normalized  by 
using  Eq.  59. 
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NORMALIZED  CORRELATION  DISTANCE 


FIG.  3  THE  FUNCTIONS  Ki(x)  AND  K2(x) 


TABLE  1  SAMPLE  VALUES  OF  Ki  AND  K2 


X 

K1  (x) 

K2  (x) 

0 

2.000000 

1 .000000 

1 

1 .723824 

0.922369 

2 

1 .295799 

0.788175 

3 

0.902805 

0.650  4  54 

4 

0.585233 

0.525517 

5 

0.344280 

0.418094 

6 

0.169415 

0.328420 

7 

0.047631 

0.254962 

8 

-0  .033259 

0.195599 

9 

-0.083592 

0.148139 

10 

-0.111697 

0.1 10546 

11 

-0.1  24076 

0.081028 

12 

-0.125690 

0.058053 

13 

-0.120  249 

0.040338 

14 

-0 .1  10477 

0.026820 

1  5 

-0.098337 

0.016631 

16 

-0.085208 

0.009063 

17 

-0.0  720  3  5 

0.003545 

18 

-0.059438 

-0.000382 

19 

-0.047799 

-0.003083 

20 

-0.037332 

-0.004850 

21 

-0.028127 

-0.005914 

22 

-0.020191 

-0.006456 

23 

-0.01 3475 

-0.006618 

24 

-0.007892 

-0.006509 

25 

-0.003340 

-0.006215 

26 

0.000295 

-0.005800 

27 

0.0031 28 

-0.005314 

28 

0.005268 

-0.004792 

29 

0.006819 

-0.004260 

30 

0.007875 

-0.003739 

31 

0.008521 

-0.003240 

32 

0.008834 

-0.002773 

33 

0.008882 

-0.002343 

34 

0.008723 

-0.001951 

35 

0.008408 

-0.001600 

36 

0.007978 

-0.001 288 

37 

0.007468 

-0.001013 

38 

0.006905 

-0.000774 

39 

0.006314 

-0.000567 

40 

0.005712 

-0.000391 
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FIG  4  SPATIAL  CORRELATION  FUNCTION  OF  WIND-GENERATED  OCEAN  SURFACE  WAVES,  DERIVED  FROM  THE  NEUMANN  SPECTRUM 
The  normalized  correlation  distances  are  2g/U^  times  the  actual  distance.  ^Kl  's  the  downwind  direction,  j^j  the  crosswind  direction. 


4 . 3  The  Pierson-Moskowitz  Spectrum 


Now  the  variable  c  is  changed  bo  z  by  means  of  the  relation 


z 


B(o0/a) 


[Eq.  65] 


as  in  Sect,  3*3*  Substitution  of  Eqs.  22,  36  and  58  into  Eq,  56 

gives  expressions  similar  to  Eqs.  61  and  62: 


p>j  C  %  ^  1  ^  ^  ( *njj  )  ^2  (  ^  ]  * 

N(  pN,  cp,  0)  =  cos2cp  L1  (  pN)  +  ( sin2  cp  -  cos2cp)  L2  (  PN)  ; 

the  functions  L  and  L  are  defined  as  follows: 

\  £ 

L,(pN)  =2  J  dz  e~7  J0(ipN  si  z-5) 


M 


O  00 


dz  e 


-  z 


1  PN B2  z  2)/(i 


pnb 


1 

'2). 


[Eq.  66] 

[Eq.  67] 

[Eq.  68a] 

[Eq.  68b] 


The  results  of  the  numerical  calculation  of  and 

presented  in  Fig.  5  and  Table  2.  Knowledge  of  these 
is  sufficient  to  simplify  the  calculation  of  N  from 
Figure  6  gives  an  impression  of  the  surface  N( 
the  Pierson-Moskowitz  spectrum. 


L2  are 
functions 
Eqs.  66  or  67* 
,  0)  for 


4 . 4  Down-Wind  and  Cross-Wind  Correlation 


From  Figs.  4  and  6,  the  function  N(?^,  0)  can  only  be 

studied  qualitatively,  because  these  figures  try  to  depict  a 
three-dimensional  surface  in  a  two-dimensional  plane.  It  is 
therefore  worth  considering  two  cross-sections  of  N(^,  T^,  0), 

namely  the  down-wind  correlation  function  N(^,  0,  0)  and  the 

cross-wind  correlation  function  N(0,  0)  .  Formulae  for  these 

functions  are  readily  obtained  from  Eqs.  61  and  66: 
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NORMALIZED  CORRELATION  DISTANCE 


FIG.  5  THE  FUNCTIONS  L,(x)  AND  L2(x) 


TABLE  2  SAMPLE  VALUES  OF  L,  AND  L2 


X 

L^x) 

La  (x) 

0 

2,000000 

1.000000 

1 

1.683109 

0.910138 

2 

1  .20  7823 

0.760491 

3 

0.765871 

0.606300 

4 

0.410755 

0.466143 

3 

0.150207 

0.346600 

6 

-0.025481 

0.248841 

7 

-0.132343 

0.171491 

8 

-0.186940 

0.112086 

9 

-0,203925 

0.067769 

10 

-0.195557 

0.035761 

11 

-0.172261 

0.013510 

12 

-0.141542 

-0.001208 

13 

-0.108166 

-0.010259 

14 

-0.075907 

-0.015168 

13 

-0.046841 

-0.017151 

16 

-0,022416 

-0,017162 

17 

-0.002859 

-0.015920 

18 

0.011061 

-0.013957 
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FIG  6  SPATIAL  CORRELATION  FUNCTION  OF  WIND-GENERATED  OCEAN  SURFACE  WAVES,  DERIVED  FROM  THE  PIERSON-MOSKOWITZ  SPECTRUM 
The  normalized  correlation  distances  are  2g/U2  times  the  actual  distance.  S  |s|  is  the  downwind  direction,  the  crosswind  direction. 


a) 


Neumann  Spectrum 


N  (  ?N  >  0  ,0)  =  -K2(?n) 

N(0,  ,0)  =  K2(ttn)  . 

b )  Pierson-Moskowitz  Spectrum 

N(Sjj,  0  ,0)  =  lx(?n)  -  ls(en) 
N(0,  ITjj  ,0)  =  . 


[Eq.  69a] 
[Eq.  69b] 

[Eq.  70a] 
[Eq.  70b] 


The  corresponding  curves  and  the  difference  between  them  are 
presented  in  Figs.  7  and  8.  We  note  that  the  correlation  in  the 
cross-wind-direction  [Fig.  8]  is  better  than  in  the  down-wind 
direction  [Fig.  7]* 


4 . 5  Discussion 

The  difference  between  the  spatial  correlation  functions  obtained 
from  the  two  energy  spectra  under  study  is  difficult  to  see  from  a 
comparison  between  Figs.  4  and  6.  This  difference  is  therefore 
plotted  in  Fig.  9* 

A  more  quantitative  idea  about  the  differences  can  be  obtained 
from  Figs.  7  and  8.  The  maximum  value  turns  out  to  be  about 
0.1226,  occurring  at  %  =  5,  = 0 . 

The  fact  that  the  correlation  distance  could  be  normalized  with 
respect  to  the  wind  speed  confirms  the  conclusion  reached  for 
the  time-correlation  function,  namely  that  the  fully-developed 
sea  changes  only  its  scale,  and  not  its  shape,  when  the  wind  speed 
is  changed  and  a  new  equilibrium  is  found. 
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FIG.  7  SPATIAL  CORRELATION  FUNCTION  OF  WIND-GENERATED  OCEAN 
SURFACE  WAVES  iN  THE  DOWN-WiND  DiRECT  ON 
The  normal  zed  d  stance  equals  2g/U^  t'mes  the  ac*ual  d  s'ance 


1. 


o) 


0.5- 


Neumann  Spectrum 
Pi er son-Mo sko witz 
Dif  ference 


30  40 


Spectrum 


— 0  #  5- 


NORMALI ZED  DISTANCE 


FIG.  8  SPATIAL  CORRELATION  FUNCTiON  OF  WIND  GENERATED  OCEAN 
SURFACE  WAVES  IN  THE  CROSS  W'ND  DIRECTION 

The  normal  zed  co'rela*  an  d  stance  equals  2g/U^  t  mes  the  actual  d  stance. 
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NEUMANN  AND  THE  PI ERSON-MOSKOWITZ  SPECTRA 


CONCLUSIONS 


Wind-generated  sea  surface  waves  can  be  considered  with  good 
approximation  as  a  Gaussian  process  G(x,  t),  stationary  in  time, 
and  homogeneous  and  anisotropic  in  space.  If  the  mean  value  of 
this  process  is  set  equal  to  zero,  which  can  be  done  without  loss 
of  generality,  the  correlation  function 

h2N(p,  t)  =  E[C(x,t)  £(x+p,  t+T)] 

represents  the  process  completely  in  a  statistical  sense.  This 
correlation  function  plays  an  important  role  in  the  description 
of  an  underwater  sound  field  that  is  scattered  from  the  sea  surface. 

The  most  realistic  way  to  describe  the  surface  waves  statistically 
is  by  using  a  surface  wave  energy  spectrum  $( C,  0).  The 
correlation  function  N(p, t)  is  related  to  the  energy  spectrum 
through  a  double  Fourier  integral,  over  wave  frequency  and  wave 
direction.  The  time- correlation  function  follows  with  p  =  0, 
the  spatial-correlation  function  with  j  =  0 •  For  two  proposed 
spectral  functions  for  fully-developed  seas,  the  Neumann  and 
the  Pierson-Moskowitz  spectra,  we  have  calculated  the  correlation 
functions  N(0,t)  and  N(^,  0)  on  a  digital  computer.  It 
turned  out,  as  was  already  indicated  by  Pierson  and  Moskowitz, 
who  normalized  the  wave  frequency  and  the  fetch,  that  the  scales 
could  be  normalized  with  respect  to  the  wind  speed:  normalized 

time  equals  real  time  divided  by  wind  speed*  normalized  distance 
is  proportional  to  real  distance  divided  by  the  square  of  the 
wind  speed.  This  indicates  that  two  fully-developed  seas  of 
infinite  fetch,  with  different  constant  wind  speeds,  differ  only 
in  scale,  and  not  in  shape. 

As  the  Neumann  and  Pierson-Moskowitz  spectra  are  only  given  as 
functions  of  wave  frequency  c,  but  not  of  wave  direction  0, 
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a  relation  between  the  directional  spectrum  and  its  isotropic 
mate  had  to  be  assumed.  We  have  chosen  the  frequently  met  (cos2  9)- 
connection,  but  a  more  adequate  relation  might  be  given  by  using 
cos13  0  ,  where  p  is  frequency-dependent. 

Although  the  correlation  function  N(p,  0)  is  essentially  a 
function  in  two  dimensions,  it  could  be  expressed  in  two 
functions  that  depend  only  on  the  modulus  of  *p .  This  is  an 
important  result  for  the  description  of  the  underwater  scattered 
sound  field  when  it  comes  to  numerical  evaluation  of  its 
statistical  properties:  instead  of  having  to  represent  N(*p,  0)  as 
a  matrix  of  M  xM  elements,  two  column  vectors  of  size  M  will 
suffice.  A  simple  calculation  then  yields  N($,  0)  for  each 
spatial  point.  Sample  values  of  these  two  functions  are  included 
in  this  report,  for  M  =  40 . 

It  is  found  that  the  Neumann  and  Pierson-Moskowitz  spectra 
produce  correlation  functions  that  differ  only  quantitatively. 


33 


APPENDIX  A 


TOTAL  ENERGY  IN  A  FULIY  DEVELOPED  SEA 

The  total  energy  E  is  defined  by  Eq.  28: 

E  =  j*  da  §(  a)  • 
o 

We  shall  now  derive  F.qs.  29  and  30. 
a )  Neumann  Spectrum 

Substitution  of  Eq.  16  into  Eq.A.l  gives: 

E  =  C  |  da  o  "  exp(-2g2/a2U2 ) , 
o 

and  with 

z  =  g  JT/  aU 
this  becomes 

TT  f  U 


\S  n » 


E  =  C 


dz  z  e 


[Eq.  A . 1] 


[Eq.  A. 2] 


[Eq.  A. 3] 


[Eq.  A. 4] 


The  integral  has  the  value  -g-  J  tt  ,  which  has  been  found  in 
Ref.  9  [p.337,  no . ( 3 • 461 . 2) ] .  And  so  we  have  finally 


n  3/2  5 

E  -  (f)  ) 


[Eq.  A. 5] 


b )  Pierson-Moskowit z  Spectrum 
We  substitute  Eq.  22  into  Eq.  A.l,  and  find: 


pCO 

ag2  ,  da  exp 


-p(-er-) 

a 


-5 


[Eq.  A. 6] 
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With  the  change  of  variables 


=  e  It, 


Eq.  A,  6  reduces  to 


E=if-ai  f  dz  e-z 
P  CT04  J0 


qg‘ 


4Ba04 


or,  because  of  <j  —  g/U  , 


[Eq.  A. 7] 


[Eq.  A. 8] 


[Eq.  A. 9] 
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APPENDIX  B 


CALCULATION  OF  THE  INTEGRAL  I^g,  p) 

We  call  cp  the  angle  made  by  the  vector  "p  with  the  ^-axis* 
So  we  have : 


sin  cp=  vj  p  y  cos  cp=  ?/  p  • 
Hence  Eq.  57  can  be  written  as 


[JEq.  B  •  1] 


^  pn/2 

I T  (  a,  p)  ==  d@  cos20  cos[o  cos  (  0  -  cp)L 

-TT/2 

where  p'  =  a2  p/g  . 


[Eq.  B . 2] 


Next  we  introduce  the  new  variable  01  =  0  -  cp>  and  drop  the 

prime .  Then  we  obtain 


la(a,  p)  = 


p  tt/  2-  CP 


d0  cos2  (  0  +  cp)  cos(p'cos0)  , 


-rr/2-cp 

and  with  the  identity 


[Eq.  B.3] 


cos2  (  0+cp)  =  cos20cos2Q5+sin20sin2cp-^sin  (20)  sin(  2<p) 

*4 

1^  can  be  split  into  three  parts: 

p  tt/  2— cp 

la  (  C,T$)  =  C°s2cp  j  d0  cos2  cos(  p'  cos  0) 

—  tt/  2—  cp 


+ 


sin2qj 


P  tt/  2—  cp 

d  0  sin2 

—  tt/  2  —  cp 


0  cos  (  p'  cos  @) 


-  ^sin(2(p) 


->  tt/  2— cp 

d0sin(20)  cos(p'cos0). 

-tt/  2— cp 


[Eq.  B.4] 


[Eq.  B.5] 
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The  three  integrands  have  the  period  tt,  because  cos2  0, 

sin2  0,  sin(2  0)  and  cos(p!cos0)  have  that  period.  And  since  we 

have  to  integrate  over  an  interval  of  length  TT,  we  can  take 

the  limits  from  -tt/2  to  +tt/ 2 .  But  then  the  third  integral 

vanishes,  as  the  integrand  is  odd.  And  the  value  of  the 

two  remaining  integrals  is  not  changed  if  we  integrate 

from  0  to  tt,  instead  of  from  -tt/2  to  +tt/2.  Hence,  we 

have  : 


P 


TT 


I1  (  c*  p)  =  cos2cp  d0  cos2  0  cos(  p1  cos  0)+sin2fp  j 

°0  0 


Substitution  of  cos2  0=1-  sin2  0  gives 

r 


pTT 

.  d 


0  sin2  0  cos(  p1  cos0), 
[Eq.  B . 6] 


n  TT 


Ia  (  Cy  p)~cos2cp  j  d0  cos(  p1  cos  0)  +  (  sin2  cp  -  cos2  cp)  d  0  sin2  0cos  (  p'  cos  0)  . 

0  Lo 

[Eq.  B.7] 


The  integrals  are  readily  evaluated  with  tables  [see  Ref.  9* 
pp. 402-403*  nos. 18  and  21*]: 


r 


rr 


cos2cpJ0(  p')  +  (sin2cp-  cos: 


J,.  (  P'  ) 


[Eq.  B.8 


With  Eq.  B.l  and  p'  =  a2 p/g  into  this  formula,  Eq.  58  follows. 
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